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Magnetization Plateaux in Random Frustrated S=l/2 Heisenberg Chains 
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The S = 1/2 frustrated Heisenberg chains with bond alternation is known to exhibit a mag- 
netization plateau at half of the saturation magnetization M a accompanied by the spontanuous 
translational symmetry breakdown. The effect of randomness on the magnetization process of 
this model is investigated. First we consider the mixture of the two kinds of chains both of 
which possess the M s /2-plateau in the common interval of the magnetization field. The plateau 
at M s / 2 is found to vanish immediately if the randomness is switched on in agreement with 
Totsuka's prediction. The small plateau also appears near the saturation field due to the local- 
ization of inverted spins around the minority bond. On the other hand, if the stronger bond 
is replaced by the ferromagnetic bonds randomly, the randomness induced fractional plateau 
appears as in the nonfrustrated case. The plateau at M s /2 does not simply vanish but shifts 
and splits into two smaller plateaux. The magnetization on this plateau varies nonlinearly with 
1 — p. The physical origin of this behavior is explained based on the strong coupling picture. 
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1. Introduction 

In the recent studies of quantum many body problem, 
magnetization plateaux in one dimensional Heisenberg 
chains are attracting much attention as the field induced 
spin gap states. On the plateau, some spins are partly 
quenched by the magnetic field while the remaining spins 
form a spin gap state. 7 ' 1,2 ' 4 ) 

Oshikawa, Yamanaka and Affleck 5 ' proposed the nec- 
essary condition for the magnetization plateaus as 

Q G s(S~M) = Z = mtegeT (1) 

where Qgs is the spatial periodicity of the magnetic 
ground state, 5 is the magnitude of the spin, and M 
is the magnetization per site. It should be noted that 
Qgs is not always equal to the spatial periodicity of 
the Hamiltonian QHam- If the translational symmetry 
is spontanuously broken, Qgs = zQnam where z is an 
integer. 

Recently, the effect of randomness on the mag- 
, netization plateaux of one-dimensional quantum spin 
systems has been studied theoretically. 6 ~ 8 ) Among 
, them, the present author found the randomness induced 
plateau in the mixture of 5 = 1/2 antiferromagnetic- 
, antifcrromagnetic alternating Heisenberg chain and 
5 = 1/2 fcrromagnetic-antiferromagnctic alternat- 
ing Heisenberg chain 8 ) which models the compound 
(CH 3 ) 2 CHNH 3 Cu(Cl a; Bri_ 2; )3 studied by Manaka and 
coworkers. 9 ) Similar observation has been made by 
Cabra and coworkers 6 ) for random 5-merized chains. 

On the other hand, Totsuka 7 ) discussed the effect of 
randomness on the magnetization plateaux which are al- 
ready present in periodic systems. He concluded that 
the plateau induced by the imposed periodicity (Qgs = 
QHam) is stable against weak randomness whereas the 
plateau induced by the spontanuous translational sym- 
metry breakdown (STSB) (Qgs 7^ QHam) is unstable 
against randomness due to Imry-Ma effect. 10 ) 

Therefore we know that randomness induces the mag- 



netization plateau on one hand and on the other hand 
it destroy the plateau. It is the purpose of the present 
work to investigate the interplay of these two apparently 
contradicting aspects of randomness effect on the mag- 
netization plateau phenomenon using the density matrix 
renormalization group (DMRG) method. 11 ) 

This paper is organized as follows. In the next section, 
the model Hamiltonian is defined. Frustrated random 
chains with bond alternation are investigated in section 
3.1 and frustrated chains with bond alternation and ran- 
dom sign strong bonds are studied in section 3.2 using 
the DMRG method. The physical interpretation of the 
numerical results are also given. The final section is de- 
voted to a summary and discussion. 

2. Model Hamiltonian 

As a simplest model which exhibits a plateau accom- 
panied by STSB, we consider the 5=1/2 Heisenberg 
chain with next-nearest-neighbour interaction depicted 
in Fig. 1. The Hamiltonian is given by, 
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Fig. 1. Heisenberg chain with next-nearest-ncighbour interaction. 
The nearest neighbour interaction is spatially modulated. 



N-l N-2 

H=J2 J i S iSi+r + JSS i S i+i ( 2 ) 

i=l i=l 

where the bond alternation and randomness is intro- 
duced in the distribution of Jj. In the present work, 
we consider the following two types of randomness: 
(1) Random bond alternation 
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Fig. 2. The magnetization curves of the pure frustrated random 
Heisenberg chains with S = 0.2, (a) a = 0.5 and (b) a = 0.6. 
The system size is N = 82. 
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In the following, we set J = 1 to define the energy unit. 

In the pure case (a± = 012 = a), this model has a 
plateau at the magnetization M = M s /2 for appropri- 
ate values of a and <5 ? ' 2 ' 4 ) where M s is the saturation 
magnetization per site (M s = 1/2 in the present model). 
To explain this plateau, we have to set Qqs — 2Qnam in 
Eq. (1). Therefore this plateau is accompanied by STSB. 
Figures 2(a) and (b) show typical examples of such cases 
with 8 = 0.2 and (a) a — 0.5 and (b) a = 0.6 calculated 
by the DMRG method. 



Numerical Results 



3.1 



Case 1: Frustrated Heisenberg chain with random 
bond alternation 
From Fig. 2, it is clearly seen that the magnetic field 
intervals of the plateau region of cases (a) and (b) have 
finite overlap. In this section, we investigate what hap- 
pens if these two chains are mixed. The result is shown 
in Fig. 3. Even in the magnetic field interval in which 
both (a) and (b) have plateau, no plateau appears in the 
mixed chain. This verifies the fragility of this kind of 
plateau due to Imry-Ma effect 10 ^ as predicted by Tot- 
suka. 7 ) 

It should be also remarked that the plateau with M — 
(spin gap) remains stable. This plateau is due to the 
imposed bond alternation (Qgs — Qitam) so that this 
result is also in accord with Totsuka's prediction. 7 ) 

The structure near the saturation is also induced by 
randomness. It should be remarked that the satura- 
tion field i/gat in the mixed chain is larger than those 
of both pure chains. The saturation field is determined 
as the energy difference between the fully polarized state 
(Mtot = N/2) and the state with a single inverted spin 
(M to t = N/2 — 1). Here M tot denotes the total magne- 
tization. In the presence of the impurity bond, the in- 
verted spin can be localized around the impurity bond. 
In Figure 4, the saturation fields of the chains with a sin- 
gle impurity bond with strength Ji mp in the host chain 
with bond alternation a = a no st is presented for N — 82 
chain. The filled circles represent the saturation field of 
the chain with a no st = 0.5 and 6 — 0.2 in which a sin- 
gle strong bond is replaced by the impurity bond with 
Ji mp = 1.6 corresponding to the strength of the strong 
bond of the chain with a = 0.6. The open circles rep- 
resent the saturation field of the chain with ahost = 0.6 
and 5 = 0.2 in which a single weak bond is replaced by 
the impurity bond with J; mp = 0.5 corresponding to the 
strength of the weak bond of the chain with a — 0.5. 
The saturation field of the pure chain is easily calculated 
analytically as 
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Fig. 3. The magnetization curves of the random frustrated Hciscn- 
berg chains with 8 = 0.2, ai = 0.5 and «2 = 0.6 for various 
values of p. The system size is N = 82 and average is taken over 
40 samples. 




Fig. 4. The saturation fields H Bll t of the frustrated Heisenberg 
chains with an impurity bond with N = 82. The filled (open) cir- 
cles represent the saturation field of the chain with a = 0.5(0.6) 
with a single strong (weak) bond replaced by the bond with 
<Amp = 1-6 (0.5). The solid and dotted lines are the saturation 
fields of the pure chains with a = 0.5 and 0.6, respectively. 
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as plotted in Fig. 4 by the solid and dotted lines for 
a = 0.5 and 0.6 with 8 — 0.2, respectively. It is seen 
that the saturation field for the chain with an impurity 
is substantially larger than the pure chains for small 5 
where the localized states are formed. Corresponding 
probability density of the inverted spin | x(i) | 2 in the 
ground state with M tot = N/2 — 1 is shown in Fig. 5 for 
(5 = 0.2. This explicitly shows that the inverted spin is 
localized around the impurity bond. 

For small but finite density of impurity bonds, the 
magnetization of the most part of the chain is saturated 
at the saturation field of the pure chain. However, the 
portion of the chain proportional to the density of the 
impurity bond remains unpolarized until the saturation 
field of the chain with an impurity is reached. Therefore 
there appears a small plateau at M = M s p for small p 
and at M = M s (l —p) for small 1 — p as observed in Fig. 
3. This can be also regarded as a randomness induced 
plateau. 

3.2 Case 2 : Frustrated chains with bond alternation 
and random sign strong bonds 
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Fig. 5. The probability density | x(i) | 2 of the inverted spin in the 
ground state with Mtot = N/2 — 1 of the frustrated Heisenberg 
chains with an impurity bond. The chain length is N = 82. 
The filled (open) circles represent | x(i) | 2 of the chain with 
a = 0.5(0.6) with a single strong (weak) bond replaced by the 
bond with J imp = 1.6 (0.5) at i = 41(40) with <5 = 0.2. 



Fig. 6. The magnetization curves of the frustrated Heisenberg 
chains with bond alternation and random sign strong bonds with 
Ja = 1.5, Jp = —1.5, J = 1.0 and 5 = 0.2 for various values of p. 
The system size is TV = 82 and average is taken over 40 samples. 
The small low field magnetization for p = 1 is dur to the finite 
size effect. 

In the absence of frustration (5 = 0, this model shows 
the randomness induced plateau at M — M s (l — p) as 
pointed out by the present author. 8 ) The same phe- 
nomenon takes place in the presence of finite 8. 

The fate of the plateau at M — M s /2 is different from 
that of case 1. The plateau no more appears at M = 
M s /2 but it does not simply vanish. The plateau splits 
into two plateaux with magnetizations M pl and M p2 as 
depicted in Fig. 6. Actually, the lower one (M p2 ) of these 
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two plateaux is only visible for relatively large values of 
1 — p. Furthermore, the magnetization on both plateaux 
depends nonlinearly on 1 — p as shown in Fig. 7. 

The nonlinear dependence of M p i and M P 2 on p can 
be understood based on the strong coupling picture for 
the plateau state proposed by Totsuka 4, 7 ) in the follow- 
ing way. On the plateau, all the spin pairs connected by 
ferromagnetic bonds are polarized. This gives the contri- 
bution M s (l — p) to the magnetization which is the value 
for the randomness induced plateau. By these bonds, 
the chain is cut into segments of finite length. Consider- 
ing that the neighbouring triplet-triplet or singlet-singlet 
pairs costs energy, 4 ' 7 ) the segments consisting of 21 + 1 
strongly coupled antiferromagnetic dimers have stably 
magnetized state with I triplets pairs. The average num- 
ber of such segment is (N/2)(l — p) 2 p 2l+1 . So that these 
segments contribute to the magnetization by an amount 



N 



21+1 



(6) 



On the other hand, for the segments consisting of 21 
strongly coupled dimers, it is enevitable to have at least 
one neighbouring triplet-triplet or singlet-singlet pairs. 
Therefore these segments have stably magnetized state 
with / — 1 or I triplets pairs. The average number of such 
segment is (N/2)(l — p) 2 p 21 . Therefore they contribute 
to the magnetization by amounts 



N 

1=1 
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or 
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(7) 



(8) 



These summations are easily carried out and the stable 
plateaux are expected at 



0.4 



0.2 



-rr7^Z__..o ° 



0.5 



1-p 



Fig. 7. The 1 — p-dcpcndcnce of the magnetization on the plateau 
which reduces to the M = M B /2 plateau in the pure p = 1 limit 
for the frustrated Heiscnberg chains with bond alternation and 
random sign strong bonds. The filled and open circles repre- 
sent M p i and M p 2, respectively. The solid and dotted cureves 
represent the formulae (9) and (10), respectively. 




Fig. 8. The magnetization curves of the random g-merized Heisen- 
berg chains with q = 4, a = 0.5, J' = 0.2 and J = 1.0 for various 
values of p. The system size is N = 100 and average is taken 
over 20 samples. 




M pl = M s (l-p+ T ^— ) 



and 



M p2 = M s (l-p + 
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(9) 



(10) 



These are also plotted in Fig. 7 by solid and dotted 
curves. The coincidence between the numerical results 
and above formulae is fairly good. The slight deviation 
is attributed to the finite size effect. 

The nonlinear p-dependence of the plateau magnetiza- 
tion is in contrast to the case of the plateau in random 
g-merized chain in which 5 = and Ji are randomly 
distributed according to the distribution. 



J(l + a.i) : probability 1 — p 
J' : probability p 



(11) 



where on = a, (—a) if i = qn, (i =/= qn). The XX model 
with this exchange distribution is investigated by Cabra 
and coworkers. 6 ) Here we present the magnetization 
curve of the corresponding antiferromagnetic Heisenbcrg 



0.5 



Fig. 9. The p-dependence of the magnetization on the plateau ran- 
dom g-merized Heisenbcrg chains (open circles) with q = 4. The 
solid line is Eq. (12). 
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model using the DMRG method in Fig. 8 for q = 4. 
In this case also, the plateau at M = M s /2 in the pure 
chain shifts due to randomness. However, this shift is 
proportional to p as, 



(12) 



as predicted by Cabra et al. 6 ^ and shown in Fig. 9 by 
the dotted curve. 

As discussed by Totsuka, 4 ' the plateau due to STSB 
at M = Ms/2 in the model (2) without randomness is 
formed by the many body effect. This plateau state can 
be regarded as a kind of Mott insulator of triplet pairs. 
This feature is also reflected in the nonlinear dependence 
of the plateau magnetization on p. On the other hand, 
the plateau state at M = M s /2 of model (11) with p = 
can be understood by the local picture. Therefore the 
plateau magnetization just shifts linearly with the num- 
ber of the impurity bonds. 

4. Summary and Discussion 

The magnetization plateaux in random quantum spin 
chains are investigated by the DMRG method. In the 
random mixture of two frustrated and dimerized chains 
with different values of dimerization parameter, it is veri- 
fied that the plateau accompanied by STSB is destroyed 
by randomness even in the interval where the uniform 
chains have a plateau in common. It is also verified that 
the plateau due to imposed spatial periodicity is robust. 
These are consistent with the analytical prediction by 
Totsuka. 7 ) The small randomness induced plateau is also 
predicted near the saturation field. This plateau is due 
to the localization of inverted spins. 

For the frustrated chains with bond alternation and 



random sign strong bonds, it is found that the nontrivial 
plateau due to STSB splits into two smaller plateaux. 
The physical picture of these plateau states are clarified. 
In contrast to the plateaux without STSB, the magneti- 
zation on the plateaux varies nonlinearly with p. 

So far, the experimental study of the interplay between 
the randomness, frustration and magnetization plateau 
phenomenon has not yet been carried out. Consider- 
ing the variety of phenomena predicted analytically and 
numerically and recent progress of high magnetic field 
technique, fruitlful physics is expected in this field in the 
near future. 
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